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We study the linear diamagnetic response of a superconducting cylinder 
coated by a normal-metal layer due to the proximity effect using the clean 
limit quasiclassical Eilenberger equations. We compare the results for the 
susceptibility with those for a planar geometry. Interestingly, for R ~ d the 
cylinder exhibits a stronger overscreening of the magnetic field, i.e., at the 
interface to the superconductor it can be less than (-1/2) of the applied field. 
Even for R^$> d, the diamagnetism can be increased as compared to the pla- 
nar case, viz. the magnetic susceptibility 4irx becomes smaller than —3/4. 
This behaviour can be explained by an intriguing spatial oscillation of the 
magnetic field in the normal layer. 



1. Introduction 

A normal metal in electronic contact to a superconductor acquires su- 
perconducting properties. This phenomenon is called proximity 
One of these properties is the diamagnetic screening of an applied magnetic 
field, which has been studied in a series of experimenta l^""^" and theo- 
retical work J^' l 1 ' I '^R 

First predictions on the induced screening properties were made on the 
basis of the Landau-Ginzburg theory by de Gennes and coworkers^'. The 
proximity effect is governed by the thermal coherence length £t, given by 
y/D/T in the diffusive regime and v-p/T in the clean limit. Here v-p denotes 
the Fermi velocity and D the diffusion coefficient inside the normal metal. 
Throughout the paper, we will set h = k& = 1. The Ginzburg-Landau 
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approach is valid outside the mesoscopic regime in which £t d, where d 
is the thickness of the normal metal film. In this approach a region of width 
~ £t of the normal metal screens the magnetic field and consequently the 
susceptibility of the normal metal 4ttx = AttM/H = —^/d increases with 
the coherence length. Strictly speaking the approach of de Gennes et al. 
requires in addition a local constitutive relation j(x) ~ A(x^^, which is only 
the case under special circumstances^. The first experimental studies were 
in agreement with these early predictions^^ and the predicted screening 
was observed. However, in this experiment the mesoscopic regime was not 
reached. According to the prediction the susceptibility approaches almost 
the ideal value 4ttx = — 1. 

In 1980 Zaikin investigated the screening properties of clean normal 
metals in the clean limit. In that case the constitutive relation is complete 
non-local, viz. j ~ A(x)dx. This has the interesting consequence that the 
susceptibility saturates at low temperatures at 75% of the ideal diamagnetic 
value. The magnetic field inside the normal metal decays linearly at low 
temperatures and even changes sign to reach a maximal opposite field of 
—H/2 at the interface to the superconductor. Note, that a similar effect 
occurs in clean type-I superconductors^. 

Around the same time a series of experiments by Mota and coworkers 
found an interesting low-temperature anomaly in the magnetic response of 
cylindrical structures I ^ I SMf l 1 ? ! 20 !. At very low temperatures the diamagnetic 
signal decreased again and in the end became even paramagnetic. This so- 
called reentrant effect of the magnetic susceptibility has triggered a number 
of theoretical explanation s! 2 1 l 22 l 23 l 24 l ; however a final experimental verification 
of one proposal is still missing. 

The effect of elastic impurity scattering was the subject of a number of 
works. A finite elastic mean free path reduces the range of the current-field 
relation and changes the screening properties drastically^. Remarkably, 
impurities can enhance the screening ability at some temperatures and even 
have an effect if the mean free path exceeds the thickness of the normal layer. 
This quasiclassical description^! was shown to agree with the experimental 
data for temperatures above the reentrance regime®. In a number of other 
works the nonlinear magnetic properties have been discussed numerically^ 
and analytically^. Other works addressed the effect of a non-ideal interface 

In this work we investigate the effect of a cylindrical geometry on the 
magnetic screening properties. A sketch of the system is shown in the left 
panel of Fig. ^ We will consider a superconducting cylinder covered by a 
normal-metal layer. The superconducting core has radius R; the normal- 
metal layer has thickness d, i.e., the total radius of the cylindrical sample is 
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Fig. 1. Left panel: a superconducting cylinder surrounded by a normal-metal 
layer. A magnetic field is applied in parallel to the axis of the cylinder. Right 
panel: parameterization of the averaging procedure over the Fermi surface. 

R + d. Both the superconductor and the normal metal are assumed to be in 
the clean limit. A magnetic field is applied parallel to the axis of the cylinder, 
and the task is to calculate the magnetic susceptibility of the normal-metal 
layer (we assume an ideal diamagnetism of the superconductor). 

We assume the pair potential A to have a step-like dependence on the 
radial coordinate, 

A(r) = A@(R-r) , (1) 

and we will neglect self-consistency. This is justified since the superconductor 
is assumed to be much thicker than the coherence length and the magnetic 
field is assumed to be much less than the critical field. A small suppression 
of A at the interface would not lead to qualitative changes. 

Throughout this paper, we will use the symmetric gauge for the vector 
potential 

A(r) = A(r)e v , H(r) = He z Q(r - R - d) . (2) 

Here is the unit vector in direction of (p and e z is the unit vector in 
z-direction. 

The rest of the paper is organized as follows: in Section EJ we introduce 
the formalism of quasiclassical Green's functions and discuss how to solve the 
Eilenberger equation. Its solution leads to an expression for the supercurrent 
density in the normal layer, that is combined with Maxwell's equation in 
Section ED which results in an integro-differential equation for the vector 
potential. We solve this equation numerically in Section ED and calculate the 
magnetic susceptibility of the normal layer. In the planar limit, the results 
agree with earlier work, whereas in the cylindrical case we find new and 
interesting screening behavior and a non-trivial oscillating field distribution. 
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2. Solution of the Quasiclassical Equations 

The goal is to solve the clean-limit Eilenberger equation I ? 1 ! ? 2 ! 

-v F dg UJ {r, v F ) = [(u + iev F A(r))f 3 + A(r) , g u (r, v F ) } (3) 

in the cylindrical geometry described above. Here uj = irT(2n + 1) denotes a 
Matsubara frequency and A = Afi is the off-diagonal pair potential of the 
superconductor. We will assume the following boundary conditions: a fully 
transparent boundary at the NS interface, 

g 3 jR-0,v F )=gZ(R + 0,v F ) (4) 

and specular reflection at the outer boundary 

ffi(R + d,v x ,v\\)=gfc(R + d,-v ± ,v\\) . (5) 

Deep in the superconductor, we have 

where f2 = \/A 2 + u 2 . These equations have to be solved along classical tra- 
jectories, along which they become effectively one-dimensional. The general 
solution is a superposition of the homogeneous solution g^ and the rising 
resp. decaying solutions g+ and In the normal layer, R < r < R + d we 
make the Ansatz 

ij(r — R) —r](r — R) 

g{r) n = a g% + a + gle +a-g!e , (7) 

where 

r,( x ) = x + 2ie^- T A(R + x') dx' (8) 
\v Fx \ \v Fx \ Jo 

is the complex phase accumulated along the trajectory and 

1 

gh=r 3 , g± = -(nTi sign(v Fx )T 2 ) . (9) 



In the superconductor, the corresponding Ansatz reads 

2Q(r-R) 2Q(r-R) 

g(r) s = Po9 s h + P + 9 s + e ^ +P-fc.e M , (10) 



where 



9h = 7^( wr 3 + Ar i) , 9± = ^( Ar 3 - WT i ± in sign(v Fx )T 2 ) . (11) 
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We have neglected the magnetic field in the superconductor, since the pen- 
etration depth is assumed to be much smaller than the coherence length. 
Furthermore, we will neglect the possibility that a quasiparticle trajectory 
leaves the superconducting core again. 

The coefficients oq, /3o,a + /_ and fulfill the normalization condi- 

tions 

a 2 Q + a + a- = 1 , $ + /3+/3_ = 1 . (12) 
Using the boundary conditions, we find that 

/?_ = and /?o = 1 (13) 

deep in the superconductor, whereas the boundary condition at x = R + d 
in the normal conductor leads to 

77(d) -77(d) 
a + e = a_ e . (14) 

The other two coefficients can be determined from the matching condition 
at r = R leading to the final result 

= Q + ^coth(7?(d)) 2A 
°° w + ncothMd)) ' P+ w + ncoth(»y(d)) ' 1 J 

The diagonal component of the Green's function in the normal metal (T3- 
component of <? n ) is constant along a trajectory, 

_ n + ^oth(„w) 

cj + S2coth(77(d)J 

where 

"~ Lc +ie [ A(r)dr (17) 



Jc 

and 

0= VA 2 + w 2 . (18) 

Here, L£ is the length of trajectory £ and J £ the line integral along the 
trajectory. Here we have generalized our previous formula to the more gen- 
eral cylindrical case. In the planar limit R/d — > 00, our solution agrees with 
earlier worlM 

This solution determines the screening current in the normal metal layer 

via 



j(r) = -i47reiVbr ^(wfTMs^Cwf)) = -iMeN a T ^(i> F a > . (19) 

uj>0 uj>0 
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Here, Nq denotes the density of states at the Fermi energy. We note, that 
the screening current is determined by the diagonal component of g only, 
which for a given trajectory is constant along that particular trajectory. 
However, the angular average may induce a space dependence, since certain 
trajectories, which do not hit the superconductor have to be excluded from 
the angular average. This will be important for the cylindrical geometry 
considered below. 



3. Solution of Maxwell's equation 

We will solve Maxwell's equation to obtain the spatial dependence of 
the magnetic field in the normal metal from the screening-current density. 
In the planar limit R/d — > do, it is possible to give an analytical solution, 
which we discuss briefly*'. In this case, the current density in the normal 
layer turns out to be constant as a function of the spatial coordinates, 



R/d- 



OO „2L „2L 

16eN v F TV / 2 60 / 2 dip sin 2 9 cos tp (20) 
M=0 Jo Jo 

A 2 (cosh <J> + cos 4>) sin (f> 



[u> sinh + ( cosh + cos (j) )Y 

Here, 



$ = ^ = ^L (21) 

V F Vp COS 9 



is an effective trajectory length, and 

r-R+d 

4> = 4 tan 9 cos <p e / A{r)dr (22) 
Jr 

is the Aharonov-Bohm phase along this trajectory. The current density 
depends on temperature and applied magnetic field H. More precisely, there 
is a non-local relationship between current density and vector potential that 
can be understood in the following way: consider the semiclassical trajectory 
of an electron in the normal metal starting at the NS interface with velocity 
vp. It will be specularly reflected at the outer interface and eventually hit 
the superconductor, see the right panel of Fig^ There it will be Andreev- 
reflected and return as a hole along the trajectory that the electron took 
before. Such semiclassical closed orbits correspond to quantum bound states 
and are known as Andreev levels^. In the presence of a vector potential, 
the electron picks up a phase 

e \ A(r)dr = 2e^ A(r)6r , (23) 



VFr Jr 
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which will lead to an energy shift for the Andreev leveltP 

rR+d 



evF v 



R 



A(r)dr . (24) 



Since the Andreev levels are homogeneous in the normal layer, this energy 
shift leads to a spatially constant current density. 

Solving Maxwell's equation with the boundary conditions 

A(R) = (25) 

and 



* A(r < 



-_ R+d R + d 



leads to the following z-component of the magnetic field: 

B(r) = H - <iirj(T,H)(R + d-r) for R<r<R + d (27) 

B(r) = H for r > R + d . 

The field may change its sign within the normal-metal layer if 

inj(T, H)d > H . (28) 

This "overscreening" effect was first discussed by Zaikin 9 . A numerical eval- 
uation of the current density Eq. (|20|) at all temperatures provides a solution 
of Maxwell's equations for arbitrary fields. 

For a cylindrical geometry, the current density in the normal layer de- 
pends on the radial coordinate. The (^-component of the current density 
reads 

j(r) = -ieN T^2(v Fv a }{r) . (29) 

LU>0 

The average over the Fermi surface can be parameterized in a way that 
allows to easily exclude trajectories that do not meet the superconductor, 
i.e., do not contribute to the current. Instead of performing the usual angular 
average 

dn VF r 1 dcose r 2n dip 

" ~ y_! 2 J 2^"' 

it is convenient to define p = r sin ip which leads to 



(30) 
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see the right panel of Fig.^for the geometrical interpretation of p. Obviously, 
only those trajectories contribute for which \p\ < R. Such a trajectory has 
length 

Using this parameterization, the line integral over the vector potential can 
be written as 

r r R+d o 

/ A(r)dx = 2 dr'Ajr) — — = signfl ■ (33) 

J Jr V r' 2 — p 2 

We finally obtain the following expression for the current density, Eq. (|29|) . 

fR rn/2 n ^-n 2 f) 

j(r) = AeN v F TY, / dp d9 P %— 2 (34) 



LU>0 



J 



27rry / r 2 — p 2 



where 



and 



A 2 sin cf> ( cosh <3? + cos <f> ) 
( uj sinh + Q ( cosh + cos (p ) ) 2 + uj 2 sin 2 (f> 

.-MC&g, (35) 

= 4e/ drJA(/)-=?== . (36) 

Using Maxwell's equations, and to linear order in the magnetic field, we 
find the following linear integro-differential equation for the vector potential: 

p 2 sin 2 A 2 (cosh$ + l) f R+d A(r')dr' 

X 



' \/r 2 — p 2 (uj sinh $ + O ( cosh $ + 1 ) ) Jr yV 2 - p 2 
Here, the effective length of the trajectories is given by 

3 = MV(R + d) 2 -p*-VW^7) _ 

vf\ sin 9\ 

We have introduced dimensionless units; lengths are measured in units of 
the (clean-limit) coherence length at T c 
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the current density contains a dimensionless material constant, 

K = 32 e 2 N vU 2 = — ( , (40) 

IT \\ N J 

where Xjsr = (47re 2 n e /m) _1//2 is a length scale that is defined in analogy 
with the London penetration depth, but using the electron density n e of the 
normal metal instead of the superfluid density of the superconductor. 



4. Results for the magnetic susceptibility 

We will now discuss the susceptibility in the linear-response regime for 
the planar and cylindrical case. 

The numerical solution in the planar case for K = 10 and a normal-layer 
thickness of 20 £ leads to the current densities and magnetic fields shown in 
Fig. |2J The left panel shows the current density (that is constant for the 
planar case). The right panel shows the corresponding spatially dependent 
magnetic fields. For T < 0.1T C the condition (|28|) is fulfilled, and the field 
changes sign at the NS interface. Below T = 0.05T C the current density does 
not increase any more, and we only show the curve for T = 0.01T C . 

As discussed before, the current density depends on the spatial coor- 
dinate in the cylindrical case. Figure 01 shows the current density and the 
magnetic field for K = 10, both the radius R of the superconductor and the 
normal-layer thickness have been chosen as 20 £. The local current density 
can be significantly larger than in the planar case, e.g., three times as large 
at the NS interface for T —* 0. The screening effect is weaker than in the 
planar case close to the outer interface; close to the NS-interface it is in- 
creased, and the overscreening effect may be even larger than in the planar 
case. 

The susceptibility of the normal-metal layer can be obtained from the 
solution of Maxwell's equation via 

where A e {r) is the vector potential in the absence of the normal layer. In 
the planar case, the vector potential can be calculated analytically and leads 
to 

4vr X (T,F)=-2vr^^. (42) 

We will now give explicit expressions for the susceptibility in the linear- 
response regime for the planar and cylindrical case. 
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Fig. 2. Current density (left panel) and magnetic field (right panel) for a 
planar geometry with £ = d/20, R/d — > oo, if = 10, and different tempera- 
tures. The material parameter K is defined in Eq. (I4U|) . 




Fig. 3. Current density (left panel) and magnetic field (right panel) for a 
cylindrical geometry with d = R, £ = d/20, K = 10, and different tempera- 
tures. 
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In the planar case, Eq. (|20|) can be simplified assuming e J^ +d A{r)dr -C 
1 or H <C &o/d 2 . In other words, the flux enclosed by a semiclassical tra- 
jectory is less than a flux quantum <I>o; the phase picked up by the electron 
has to be less than 2tt. For experimentally relevant normal-layer thicknesses 
of d ~ l^um these conditions are fulfilled for fields H < 20G. The current 
density (|2T)|) factorizes like 

[R+d 

j(T,H) = -ej s (T) A(r)dr, (43) 

JR 

where 

oo „z 

j s (T) = 4ireN v F T V j d9 sin 2 9 tan 6 (44) 

A 2 ( cosh f ^^ 7? N ) + 1 
\ \vp cos y / 



^^h(^) + fi(cosh(^| 7 ) + l^ ' 2 ' 

The analytical solution of Maxwell's equation, Eq. ()27|) leads to the following 
expression for the susceptibility 

, , -37rej s (T)d 3 , . 

4lr ^ r) = 3 + 4,(r)rf> ' (45) 

In general, j s can be calculated only numerically. However, for T = an 
approximate analytical expression can be obtained, 

For d S> Aiv, Eq. (|43|) leads to a susceptibility that it is independent of 
geometry, 

4vr X (r = 0) = ~ . (47) 

Field expulsion is never complete; the maximal value that the susceptibility 
can reach at any temperature is 3/4 of an ideal diamagnetP 

Figure^shows the susceptibilities for different normal-layer thicknesses. 
In thicker layers, field expulsion occurs only at lower temperatures. The 
thickest layer reaches the saturation value given in Eq. (|47|) for T — > 0, since 
it fulfills d > Aat. 

For the cylindrical case, the integro-differential equation for the vector 
potential, Eq. (|37jl. has to be solved numerically. 

Surprisingly, at low temperatures and for K 3> 1, the susceptibility 
does not go to —3/4. To analyze this parameter range in a quantitative 
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Fig. 4. Susceptibility of the normal-metal layer as a function of its thickness 
d for R/d — > oo and K = 1. 

way, we start from Eq. (|37jl. We assume d <C R and concentrate on the 
contribution of the (long) trajectories with R — p ~ d. The typical length of 
these trajectories is L ~ \J~dR. They become important at low temperatures 
T -C vp/y/dR. We will assume the zero-temperature limit and replace the 
sum over oj in Eq. (|37|) by an integral. This integral converges at w — v p/L <C 
A, and the integration is readily performed. The next step is to implement 
the condition d <C R. To this end, we do the following substitutions: 

p = R- fid, (48) 

and 

r = R + zd; r' = R + z'd , (49) 

where p, z, z' are dimensionless variables with p > and < z, z 7 < 1. For 
instance, this gives \J r'2 — p 2 — > \JldR\J z' + p. Here and in the following, 
we will write A(z) instead of A(R + zd). After this substitution, Eq. 
takes the following form: 

m 2 d a A_ K fd\ l ' 2 /- d /- 1 ^ A^) 

\d) dz 2 12y/2 \RJ Jo /i Jo Z {yfp + T - x /p)^/z~Tpy/z T Tp ' 

(50) 

There is a subtle point here: the integral over p diverges at big p. Therefore, 
the right-hand side (r.h.s.) cannot describe the constant part of the current 
density correctly, however, it gives the correct inhomogeneous part. To get 
a physical result we perform the following renormalization procedure: we 
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introduce a yet unknown constant C on the r.h.s. to reproduce the constant 
part of the density and rewrite Eq. (|5U|) as 

A 2 d?A(z) _ K ( d\ 1/2 [°° J f 1 , , 

^( (vmTt-^W^ - 4 '^ 1 -^') ■ (51) 

The //-dependent expression has been subtracted from the second term on 
the r.h.s. to assure that after subtraction the term gives zero if integrated 
over z, z' . The integral over \i converges now. (We could use any other form 
of this subtraction that provides convergence of the integral since the change 
would be incorporated in C). 

We do not have to know the coefficient that relates C and A provided 
that this coefficient is sufficiently big, which is true for A <C d. In this limit, 
C has to be determined from the extra condition f dzA(z) = 0. 

The remaining task is to numerically solve Eq. (|51|) under the conditions 

f dzA(z) = ; (52) 
J o 

A A (z)\ z=l =Hd; (53) 
dz 

A(0) = . (54) 

The susceptibility is determined from this solution as 

4ir X = A(l)/Bd-1 . (55) 

Looking at Eq. (|51|) we see that the susceptibility is a function of the pa- 
rameter 

7 es l/^K(d/0 2 (d/RY/ 2 = ^f A {\/d)(R/d) 1 '* . (56) 

If 7 S> 1, 47rx goes to the known value of —3/4 and the spatial profile of A 
is determined by the magnetic energy. If 7 <C 1, the inhomogeneous terms 
dominate and Att\ goes to —1. The crossover takes place at 7 — 1. 

Figure|S]shows the susceptibility as a function of the parameter 7 defined 
in Eq. (|56[) . The universal behavior predicted by Eq. I|51|) and following 
(solid line) is compared to a full numerical solution of the integro-differential 
equation Eq. (|3*7|) for R/d = 10 8 (symbols). The two curves agree perfectly 
for 7^5; the deviation for larger values of 7 is explained by the fact that 
the condition A <C d assumed in the derivation of the universal curve is not 




Fig. 5. Susceptibility as a function of the parameter 7 denned in Eq. (|5H|) . 
Solid line: universal behavior predicted by Eqs. (|5 1 j) - (|55[). Symbols: solu- 
tion of Eqs. ((S3), (ESI, ©, and (HU for i?/d = 10 8 . 



fulfilled. The susceptibility shows a very interesting behavior. For small 
values of 7 the susceptibility deviates substantially from the planar case. In 
particular, its absolute value increases, i.e., the screening is enhanced by the 
cylindrical geometry. This is in contrast to the previously discussed behavior 
for Rk, d. The difference can be understood by looking at the magnetic field 
inside the normal metal, see Fig. H3 The spatially dependent integral kernel 
allows the existence of higher Fourier components of the current, instead of 
only the lowest-order component j =const. as in the planar case. This leads 
to characteristic spatial oscillations of the magnetic field, which can show 
several sign changes inside the normal-metal layer. In particular, the field at 
the interior interface to the superconductor is also oscillating as a function of 
the parameter 7. It is interesting to note that a similar phenomenon occurs 
in the so-called Pippard superconductors^, and has been experimentally 
observed^. 

Finally, we would like to point out the qualitative similarity of the effect 
discussed and the reduced diamagnetic response found in Ref.^^for a planar 
geometry and anomalously small impurity concentrations. In both cases, the 
origin of the effect is the absence of the contribution of very long trajectories. 
While in the cylindrical geometry the typical length scale of these trajectories 
is V dR (3> d for d -C R), the presence of impurities limits the trajectory 
length to the elastic mean free path. 
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Fig. 6. Magnetic field distribution in the normal layer for different values 
of 7. 



5. Conclusion 

We have investigated the induced magnetic screening properties of cylin- 
drical normal metal-superconductor heterostructures. The ballistic screen- 
ing leads to a peculiar non-local current-vector potential relation inside the 
normal metal, governed by the parameter 7 defined in Eq. (|56|) . For 7 > 1 
the screening is like in the planar case, viz. the current density is constant, 
the magnetic field decays linearly down to —H/2 and the susceptibility is 
4-7TX = —3/4. In contrast, for 7 < 1, the field has an oscillatory spatial de- 
pendence, including several sign changes inside the normal metal. Overall, 
the magnetic susceptibility is enhanced and reaches the ideal diamagnetic 
value —1 in the limit 7 <^ 1. 

This work was financially supported the SFB 513 Nanostructures at Sur- 
faces and Interfaces of the DFG, the Swiss NSF, and the NCCR Nanoscience. 
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